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An approach is developed for constructing simple analytical formulae accurately approximating 
solutions to eigenvalue problems of quantum mechanics. This approach is based on self-similar ap- 
proximation theory. In order to derive interpolation formulae valid in the whole range of parameters 
of considered physical quantities, the self-similar renormalization procedure is complimented here 
■ by boundary conditions which define control functions guaranteeing correct asymptotic behaviour 

| in the vicinity of boundary points. To emphasize the generality of the approach, it is illustrated by 

different problems that are typical for quantum mechanics, such as anharmonic oscillators, double- 
well potentials, and quasiresonance models with quasistationary states. In addition, the nonlinear 
Schrodinger equation is considered, for which both eigenvalues and wave functions are constructed. 
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A standard problem in quantum mechanics is how to solve approximately stationary Schrodinger equations that do 
Q^ , not possess exact solutions. In the cases of asymptotically small and asymptotically large coupling parameters, one 
■ may employ perturbation theory presenting solutions as power series in powers of the corresponding small parameter. 
However, such series are practically always only asymptotic and quickly diverge for sufficiently small expansion 
parameters. Moreover, physical quantities of interest usually correspond neither to weak-coupling nor to strong- 
coupling limits, but to an intermediate region of a coupling parameter. Thus, the problem arises of how to construct 
an interpolation formula, valid in the whole region of physical variables, when only asymptotic expansions near 
g ' boundaries are known. 

, The most known method of deriving interpolation formulae is the two-point Pade approximation [1-3]. In some 
cases the latter yields quite reasonable results. Nevertheless, the usage of this method has not become widespread 
| because of the following difficulties: 

First of all, to reach sufficient accuracy by employing Pade approximants, one needs to have tens of terms in 
perturbative expansions [1-3]. But the standard situation in physically interesting problems is when one has in hands 
only a few terms. In such a case, for the same problem one may construct different two-point Pade approximants, all 
having correct left-side and right-side limits, but differing from each other in the intermediate region by 1000% [4]. 
This clearly shows that in the case of short series the two-point Pade approximants cannot provide even qualitative 
description. 

Second, two-point Pade approximants can treat at infinity only rational powers [1-3] and are not able to describe 
other types of behaviour, for example, power laws with irrational powers or exponential functions. However, more 
complicated than the rational-power behaviour often occurs in physical problems. For instance, exponential behaviour 
at infinity is constantly exhibited by wave functions. 

Third, interpolation between two different expansions, by using two-point Pade approximants, can be accomplished 
solely, when these two expansions have compatible variables [1-3]. For example, even for a such a simple problem as 
the anharmonic oscillator, the eigenvalues in the weak-coupling and strong-coupling expansions have incompatible 
variables [5]. 

Fourth, there exists the well known and annoying problem of appearance of poles in Pade approximants, which 
results in unphysical singularities [1-3]. Eliminating such singularities in two-point Pade approximants is often 
impossible because of restrictions that are imposed by prescribed boundary conditions. 

Finally, Pade approximation is rather a numerical technique, but we keep in mind analytical approach that would 
combine relatively simple representation for physical quantities with their good accuracy. The advantage of having 
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analytical expressions, as compared to just numbers obtained from a numerical procedure, is in the convenience of 
analyzing such expressions with respect to physical parameters entering into them. 

In the present paper we develop an analytical approach for deriving interpolation formulas, which is free of the above 
deficiencies of Pade approximation. This approach works well when just a few terms of asymptotic expansions arc 
available; it successfully sews power-law with exponential asymptotic behaviour; it does not have at all the problem 
of compatibility; no unphysical poles arise; it combines analytical representation with good accuracy. We illustrate 
the approach by several quantum-mechanical problems that are usually considered as typical touchstones for any new 
method. These problems include calculation of energy levels for different anharmonic oscillators, for the Hamiltonians 
with double-well potentials, and for quasiresonance models. Each of these problems has its own specific calculational 
difficulties (for review see Refs.[5-7]). This is why it is important to show that all of them can be treated by the same 
approach. Moreover, we demonstrate that the same method is applicable to the nonlinear Schrodinger equation, for 
which we find both the energy and wave function of the ground state. The latter example is interesting not only as 
an illustration of wide possibilities of the approach, but it is important for practical purpose, being related to the 
description of Bose-condensed particles in traps. We carefully compare the properties of the wave function we have 
derived with those of the Thomas-Fermi and variational-Gaussian approximations. The analysis proves that our wave 
function provides the best approximation for a solution to the nonlinear Schrodinger equation considered. Analytical 
expressions for the wave function of a vertex filament are also constructed, being in good agreement with numerical 
data. 



II. SELF-SIMILAR INTERPOLATION 



Assume that we are interested in finding a function f(x) in the interval x\ < x < X2- The latter can be finite 
or infinite. Let equations defining the function f(x) be rigorously unsolvable, so that only perturbative asymptotic 
expansions near the boundaries can be derived: near the left boundary, 

f(x) ~ p k (x,xi), (1) 

and near the right boundary, 

f(x) ~ pk{x,x 2 ), x^x 2 -0, (2) 

where k = 0, 1, 2, ... For the time being, we do not specify the physical nature of the function f(x) and its variable x, 
since the general scheme does not depend on these specifications. 

In this section we develop such a general scheme for constructing approximations to the function f(x), so that these 
approximations, interpolating between the asymptotic expansions (1) and (2), could be valid in the whole region x± 
< x < X2- The approach we develop is based on the self-similar approximation theory [8-13] in its algebraically 
invariant formulation [14-16]. Here we show how to construct self-similar approximations so that they be compatible 
with the asymptotic boundary conditions (1) and (2). Since all theoretical foundation and basic technical details of 
the method have been expounded in our previous papers [8-16], we do not repeat them here but only delineate the 
scheme of the approach adapting it to the considered problem of interpolation. 

Let us take an asymptotic expansion, like (1) or (2), in the vicinity of a point x i7 with i = 1,2. Define the algebraic 
transform 

P k (x,s,Xi) = x s p k (x,Xi), (3) 

whose inverse, evidently, is 

p k (x, Xi) = x~ s P k (x, s, Xi). (4) 
Introduce an expansion function x((p,s,Xi) by means of the equation 

Po(x, S,Xi) — ip, X = x(ip,8,Xi). (5) 

Substituting this expansion function into Eq. (3), we obtain 

y k (tp,s,Xi) = P k {x(<p,s,Xi),s,Xi). (6) 

Transformation inverse to (6) reads 
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Pkix.s.Xi) = y k (P Q (x,s,Xi),s,Xi). (7) 

The family {yu] of the endomorphisms defined in (6) is called [11-13] the approximation cascade, because its trajectory 
{yk(f, s, Xi)} is bijective to the sequence of approximations {Pk(x, s, ir,)}. The cascade velocity can be given by the 
finite difference 

v k (ip,s,Xi) = y k (<p, s,Xi) - y k -i(ip,s,Xi). (8) 
The evolution equation, written in the integral form is 



dtp 



J 

JPk-i 



V k (<P,S,Xi) 



(9) 



where P k = P k (x, s,Xi), the upper limit P£ = P k {x, s,T,Xi) is a self-similar approximation corresponding to a 
quasifixed point, and r is an effective time necessary for reaching this quasifixed point. The latter, in accordance with 
the inverse algebraic transform (4), yields 

p* k (x, S, T, Xi) = X^ S P k (x, S, T, Xi). (10) 

To illustrate these steps, consider an asymptotic expansion 

k 

Pk (x,0) = Y / a n x n (11) 

n=0 

in the vicinity of x\ = 0. Then, accomplishing the described procedure, for Eq. (10) we find 

s/fc 



sa n 



p* k (x,s,T,0) 

An important particular case is when s — ► oo, then Eq. (12) gives 



(12) 



lim p%(x,s,T,0) =p fe _i(x,0)exp ( — rx k ) . (13) 

This shows how exponential functions naturally appear in our method, together with the radical expressions of type 
(12). 

An expression p* k , given either by (12) or by (13), as is seen, is a function of a lower-order series Ph-i, 

Pk = Fk(pk-i). (14) 

Analogously to the way by which we have come from an asymptotic series p k to the renormalized expression p* k , we 
can renormalize p k -i entering into relation (14), which gives 

Pk = Fk(pl-i) = FkiFk-Apk^))- (15) 

Repeating such a renormalization k times, we come to 

P * k -* = F k (F k _ 1 (...F 1 (p ))...). (16) 

At each n-step of renormalization (14), two parameters, s„ and r„, arise in the resulting expression, according to 
(12). Therefore, the fc-times renormalized quantity (16) contains 2k such parameters, 

Pk"* = Fk(x,s k ,f k ,x t ), (17) 

where the short-hand notation 

S fe = {si,S2,-..,S k }, f k = {n,T2,...,Tk} 

is used. 
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The sets Sk and fk are to be defined so that the renormalization procedure would converge to a function satisfying 
the boundary conditions (1) and (2). Suppose that we started from a series Pk{x, Xi) written for an asymptotic region 
of Xi, with i = 1,2. Renormalizing this series 2k times, we get (17). In order that the renormalized expression (17) 
could satisfy the correct asymptotic behaviour at another boundary point Xj, with j ^ i, we have to require the 
asymptotic condition 

F£(x,Sk,fk,Xi) — > pk{x,Xj), x— >Xj. (18) 
Condition (18) defines the control sets 

Sk = s k {x), f k =fk{x) (19) 

of control functions si(x), S2{x), . . . , Sk(x), and ti(x), T2(x), . . . , Tk(x). Substituting these control functions into (17), 
we obtain the final self-similar approximant 

fk(x,Xi) = FZ{x,s k {x),fk{x),Xi). (20) 

Control functions are called so because of their role of controlling convergence of the procedure to a function having 
the desired properties [17]. In general, these functions are, really, functions of x, although in particular cases they 
can become just parameters. In the latter case, they can be called control parameters. 

In order to make the above procedure transparent, let us consider a typical case of two asymptotic expansions at 
xi = and X2 — oo. Assume that at the left boundary we have a sequence 

Pi(x, 0) = a + aix, P2(x, 0) = a + a\x + a 2 x 2 , . . . (21) 

of perturbative expansions Pk(x, 0), and at the right boundary, a sequence 

Pi(x,oo) = Ax n , p 2 (x,oo) = Ax n + Bx"\... (22) 

of asymptotic expressions Pk{x, oo), with n > m. Starting from pi(x, 0), according to (12), we get 

pl(x,s,T,0)=^ s --^x^ . (23) 

As the asymptotic boundary condition (18), we have 

p\(x, s,t,0) ->pi(x,oo), x^oo, (24) 
with pi(x, oo) given by (22). Condition (24) holds true if and only if 

/ /IN 1 /" 

a ( 



s = —n, t = n — — . (25) 

Therefore, the first-order self-similar approximant, defined by (20), becomes 

A*(x,0) = (al /n + AV n x) n . (26) 
Similarly, starting from p2(x,0) given in (21), we find the twice renormalized expression 

f 9 1 ~ S2/2 

F*(X,S 1 ,T 1 ,S 2 ,T2,0) = \pUx,Sl,Tl,0)]~ V " ~ 1+2/s ^ ' ( 27 ) 

I s 2a J 

Imposing the asymptotic boundary condition 

F 2 *(x, si,ti,S2,t 2 ,0) -»pi(ar,oo), x -» oo, (28) 

we obtain 

na ( A\ 2/n 

s 2 = -n, t 2 = — — . 29) 
2a 2 \ a o/ 
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Employing (29) for (27), we have 



f; (x, Sl , n , 82, t 2 , 0) = { \pl (x, Sl ,n, 0)] 2 /" + A 2 ' n x 2 } 



n/2 



The boundary condition 



is satisfied provided that 



where the notation 



F 2 *(o;,si,ti,S2,T2,0) -^p 2 {x,oo), 



si 



a f 2/nq 

n = -si— — 



00, 



2B\ Vfl 



<7 = 2 + m — n<2 (n > m) 
is used. With the control parameters given by (32), the function entering into (30) writes 



pl(x,si,n,0) = 



a 2/n q + A 2, m 



2B 
nA 



1/9 



-i nq/2 



(30) 



(31) 



(32) 



(33) 



(34) 



Combining (30) with (34), we obtain the second-order self-similar approximant 



,-' n * + A 2/nq ( 

ft A 



1/9 



n/2 



A 2 ' n x 2 



(35) 



defined in Eq. (20). In the same way, we may proceed farther calculating a fc-order self-similar approximant. 

To complete this calculational procedure, we need to answer the following question. Assume that we have two 
asymptotic expansions near two boundary points. We may start from one of these expansions, say Pk(x, Xi), imposing 
the boundary condition, as in (18), at another boundary point, in this case at x 2 , and thus obtaining the self- 
similar approximant f£(x,xi). The same procedure could be accomplished, starting from pk{x 1 X2) and imposing 
the boundary condition at x\, thus getting f£(x,X2)- The question that arises is which of these two approximants, 
f£(x,xi) or f£(x,X2) is expected to be more accurate? 

The answer to this question can be done from the point of view of stability analysis [11-13]. To this end, let us 
take an expansion Pk(x, Xi) near a point Xi, with i = 1,2. Suppose that Po(x, Xi) depends on x. If p®{x, Xi) does not 
depend on x, we have to take pi(x,Xi). Define the expansion function x(ip,Xi) by the equation 



Po(x,Xi 



x(ip,Xi). 



Introduce 



yk(<P,Xi) =p k (x(ip,Xi),Xi), 



(36) 



(37) 



being a trajectory point of an approximation cascade {y l k } formed by the family of endomorphisms from (37). The 
stability of the cascade trajectory is characterized by the local multipliers 



whose images in the x-space are given by the local multipliers 

5pk(x,Xi) _ dp k (x,Xi)/dx 



m k (x,Xi) = (j,k(Po(x,Xi),Xi) 



Spa(x,xi) dpo(x,Xi)/dx' 



(38) 



(39) 



The smaller absolute values \m,k(x,Xi)\ of the multipliers correspond to the more stable trajectory of the associated 
cascade, and the higher stability implies the better convergence property of the related sequence of approximations 
[12,13]. Therefore, in the asymptotic boundary condition (18), we must choose that asymptotic expansion pk{x,Xj) 
which corresponds to the more stable cascade trajectory. 
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If two multipliers, rrik(x,x\) and rrik(x,X2), are equal or close to each other, then we cannot decide apriori which 
of the self-similar approximants, fjl(x, x\) or f£(x, X2) is preferable. In such a case, it is logical to define the average 
self-similar approximation 

r k (x) = l[r k (x,x 1 )+r k (x,x 2 )]. (40) 

Usually, one of the approximations f£(x,Xi) where i = 1,2, lies below, and another above the exact function f(x). 
In such a case, the errors of these approximants compensate each other, essentially improving the accuracy of the 
average approximant (40). 



III. ANHARMONIC OSCILLATORS 



We start illustrating our interpolation approach with the models of one-dimensional anharmonic oscillators de- 
scribed by the Hamiltonian 

H = - 1 -^ + 1 -x* +9 x- (41) 

in which the space variable x e (—00, +00), the coupling parameter g e [0, 00), and the power to > 4. These models 
are classical touchstones from which everyone starts considering a new method. 

Let us be interested in finding the ground-state energy e(g) as a function of the coupling parameter g. For this 
function, the asymptotic expansions in the weak- and strong-coupling limits are known. In the weak-coupling limit, 
perturbation theory gives 

k 

e k (g,0) = J2 a "9 n (ff-0). (42) 

n=0 

This series strongly diverges for any g ^ 0, since the coefficients a n grow like n! as n — > 00 [18,19]. The coefficients 
a n are, of course, different for different types of oscillators, depending on to. However, for the sake of simplicity, we 
do not use the double indexation. In the strong coupling limit, one has [20] the expansion 

e k (g, 00) = £ ^(i-2„)/ (ro+ 2) {g _ oo). (43) 

71=0 

Here again the coefficients A n depend on to, that is, on the kind of oscillator. Not marking this dependence explicitly 
will not lead to confusion, since different kinds of oscillators will be considered separately. 

A. Quartic Oscillator 

Start with the quartic oscillator with to = 4. For the first several coefficients of the weak-coupling series (42), one 
has [18,19] 

1 3 21 333 

The coefficients of the strong-coupling expansion (43) have been computed by many authors, starting from Hioe 
and Montroll [21]. One of the most accurate computations have been accomplished by Weniger [22]. The values 
of the strong-coupling coefficients are A = 0.667986, A x = 0.143669, A 2 = 0.008628, A 3 = 0.000818, A 4 = 
0.000082, A b = 0.000008. 

Following the approach described in Sec. II, we may start from the weak-coupling series (42) and define control 
functions from the asymptotic condition (18) with e k (g, 00) given by the strong-coupling expansion (43). In the 
second order this gives 



1/6 

4(9,0) 



9ar 
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with the control parameter 



TO 



3aQai \6A Ai 



1/4 



0.66046. 



As a result, we obtain 



where 



e* 2 (g,0) = k(l + C.g) 4 / 3 + ^ 



1/6 



(44) 



C = 



6A A } 



1/4 



In the same way, starting from the strong-coupling expansion (43) and defining control functions from asymptotic 
condition (18) at the left boundary, we find 



e%(g,oo) = 



4+{l&ata\+Alg 2 ^) 112 g 



1/4 



(45) 



Numerical calculations show that both approximants (44) and (45) are close to each other. The accuracy of a self- 
similar approximant e£(<7, Xi) can be estimated by comparing the values, it gives for different coupling parameters g, 
with precise numerical calculations accomplished by Hioe and Montroll [21] for g in the interval 0.02 < g < 20000. 
The maximal percentage error of the left approximant (44) is —2.9% occurring at g = 0.3, and the largest error of 
(45) is 4.2% at g = 2. For all g, the left approximant (44) lies below the exact values of the energy, while the right 
approximant (45) is above the exact values. The average self-similar approximant 

4(g) = I [ e 2(9>°) + e 2(.9>°°)] 



has the maximal error of 1.4% at g = 2. 

As we see, quite simple analytical expression provide sufficiently good accuracy, with the maximal error around one 
percent. As far as the structure of pcrturbative series for the quartic anharmonic oscillator is analogous to that of 
series for the so-called if 4 model of quantum field theory [23] , we may hope that for the latter one also could construct 
analogous self-similar approximants. 

B. Sextic Oscillator 

The sextic oscillator (m = 6) is interesting being a borderline case between the models whose perturbative series 
are Pade summable and those whose series cannot be summed. For the sextic oscillator, Pade approximants converge 
so slowly that they are computationally useless [5,24]. 

Employing the approach of Sec. II, we use the coefficients A = 0.680703, A 1 = 0.129464, A 2 = 0.005512. A 3 = 
0.000328, Ai = 0.000018, A 5 = 0.000001. We find the first self-similar approximants, from the left, 



e* 1 (g,0) = ^(l + WA 4 g) 



1/4 



and from the right, 



e\{g^)= l -{l + AAlg^) 



1/2 



(46) 



(47) 



The maximal error of Eq. (46), with respect to accurate numerical results [5] that can be treated as exact, is about 
— 11%, and the maximal error of (47) is about 8%. 
In the second order, we find 



e* 2 (g,0)^-[(l + 2Cg)^ 2 + (2A fg 2 



1/8 



(48) 



where C = 3.428, with the maximal error —6%. The right approximant has a comparable accuracy. 
The third-order left approximant is 
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1 f r , 1 5/4 1 1/12 

el{g,Q) = -\[(l + 2B l gfl 2 +AB 2 g 2 \ + (2A ) 12 . 9 3 } ■ 

where B\ = 4.831 and B 2 — 9.352. The maximal error is around —4%. 
The fourth-order approximant can be written as 



(49) 



(l + 2C\g) 3 / 2 +4C 2 g 2 + 8C 3 g 3 



5/4 



7/6 



1/16 



16 4 



+ (2A ) ib 5 



(50) 



with Ci = 6.078, C 2 = 18.143, and C 3 = 22.322. The maximal error is about -3%. 
In the fifth order, we find 



e£(ff,0) = - 



(1 + 2D ig f' 2 + AD 2 g 2 



5/4 



-,7/6 



8^33 J 



-1 9/8 



16D 4 g 4 



1/20 



(24>) 2 V 



(51) 



where D 1 = 7.215, L> 2 = 28.848, D 3 = 56.001, and D 4 = 49.39. The maximal error is -2.5%. 
For the sixth order, we obtain 



e* 6 (9,0) 



(1 + 2K l9 f/ 2 + 4K 2 g< 



n5/4 



+ 8K 3 g< 



7/6 



+ l6K 4 g 4 



9/8 



-1 11/10 



+ 32X 5 ff 5 



+ 



+ (2^o) 2 V 



1/24 



(52) 



with the coefficients K x = 8.256, AT 2 = 41.122, K 3 = 109.122, K 4 = 153.119, and K 5 = 104.156. The maximal error 
of approximant (52) is —2%. 

Eqs. (46)-(52) show that the accuracy of the self-similar approximants improves with increasing order. To demon- 
strate that there is uniform numerical convergence for all g, we present in Table I the percentage errors e* k (g,Xi) of 
the corresponding approximants e* k (g,Xi), as compared to exact values e(g). The accuracy in each order can also be 
improved by defining the average approximants (40). We show this for the case of the approximant 



e*i(9) = ll4(9,0) + e* 1 (g,oo) 



(53) 



whose errors are also presented in Table I. 
C. Octic Oscillator 

The case of the octic oscillator (m — 8) is important to consider remembering that Pade approximants are not 
able to sum the corresponding perturbation series [24,25]. As we show below, in our approach we obtain a series of 
self- similar approximants exhibiting uniform numerical convergence. 

Here we use the following coefficients A = 0.704046, A x = 0.120626, A 2 = 0.004168, A 3 = 0.000188, A 4 = 
0.000007, A 5 = 0.000001. The first-order left approximant reads 



el(g,0) = ^(l + 32Alg) 



1/5 



(54) 



1/2 



while the right approximant is 

el{g^)= l -(l + AA 2 g 2 ^- . (55) 
Comparing this with numerical results [5], we find that the maximal error of (54) is about —13% and that of (55) is 



For the second-order left approximant we have 



e* 2 (g,0) = ll(l + 2Cg) s / 5 + (2A ) w g 2 



1/10 



(56) 
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with C — 5.944, the maximal error being —8%. 
In the third order, we get 



1 fr , n 13/10 ~] 1 / 5 

\{g, 0) = - | [(1 + 2B ig f/* + AB 2 g 2 \ + (2A ) 15 g 3 J 



(57) 



where B 1 = 8.671 and B 2 = 26.807. The maximal error is -6%. 
The fourth-order approximant is 



el (5,0) 



1 



(1 + 2C l9 f/ 5 + 4C 2 ff s 



-i 13/10 



+ 8C 3 g 3 



6/5 1 V20 

+ (2Ao) 20 9 4 



(58) 



with Ci = 11.151, C 2 = 55.077, and C 3 = 104.667. The maximal error is -4.5%. 
The fifth order approximant writes 



(1 + 2D l9 f/ b + AD 2 g 2 + 8D 3 g 



13/10 



6/5 



16D 4 g A 



23/20 



1/25 



(24>) 2 y 



(59) 



where L>i = 13.443, D 2 = 91.126, D 3 = 282.775, and D 4 = 377.013. The maximal error is -3.5%. 
For the sixth order, we obtain 



4(9,0) = -{ 



(l + 2K ig ) 8 ^ + 4K 2 g 2 + 8K 3 g 



13/10 



6/5 



-i 23/20 



+ l6K 4 g 4 



+ 2,2K 5 g- 



28/25 



+ 



-(2A ) 3 V 



1/30 



(60) 



where K x = 15.508, K 2 = 133.486, K 3 = 581.021, K 4 = 1274, and K b = 1291. The maximal error is -3%. 

As we see, in our approach there is no principal difference between the types of oscillators, whether it is quartic, 
sextic or octic; for each of them we can easily construct a uniformly convergent sequence of self-similar approximants. 
The accuracy of the latter in each order can be essentially improved by composing average approximants, as in (53). 
The errors of the obtained approximants are collected in Table II. 

Let us emphasize that our aim here was to derive analytical formulas. The approximants we have constructed are 
easier to use than more complicated expressions that follow from renormalizcd perturbation theory [17], in which 
control functions are introduced into a zero-order Hamiltonian [6,7,26-31]. This especially concerns the sextic and 
octic oscillators. We think, that the possibility to have simple and accurate formulas, valid for the whole range 
of coupling parameters is an advantage of our approach that could be a useful tool for analysing the properties of 
quantum dots [32-34]. 



IV. DOUBLE-WELL OSCILLATORS 



Models with double-well potentials are notorously known to be difficult for approximate treatment. For instance, 
perturbation theory in this case results in series that are not Pade summable. At the same time such potentials are 
quite common for various problems encountered in physics and chemistry (see discussion in Ref.[35]). 

One of the difficulties of dealing with double- well models is that the corresponding physical quantities, as functions 
of the coupling parameter, can display not two characteristic regions of behaviour, that is, the weak-coupling and 
the strong-coupling regions, but a third region, intermediate between weak and strong coupling. This behaviour is 
similar to that of some models of quantum field theory where in the transition region instanton effects are crucial, 
bridging the weak and strong coupling limits [36,37]. 

A. Zero Dimensional Model 

Let us, first, consider the so-called zero-dimensional double-well model whose free energy is written as 

f(g) = -\nZ(g), (61) 
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where 



Z(g) = ^=[ exp {-#(*)} 

V 7 *" J-oo 



plays the role of a partition function with the Hamiltonian 

H(x) = -x 2 + gx 4 , g>0. 

The latter has a maximum H(0) = at x = and two minima H(±a) = —1/4(7 at x = ±a = ±l/^/2g. 

Direct use of perturbation theory, in powers of g, to the free energy (61) is impossible, since f(g) — ► — oo as g 
Thence, a special procedure is necessary. To this end, we define the trial Hamiltonians 

H± (x) — lo 2 (x ± a) 2 — uq, 

in which 

1 1 



and w is a trial parameter to become latter a control function. The partition function (62) can be written in the 

i r+oc i r+oc 

Z(g) = -^ c X p{-H + (x)-AH + (x)}dx + -^ exp {-H_(x) - AH_(x)} dx, 

2 V 7T J-oo 2 V 7T J-oo 

where 

AH± = H(x) - H±(x). 
The free energy (61) is expanded in powers of AH±(x), with the zero-order term 

F (g,w) = lnw - u , 

the first-order term 

F l(5 ,.) = ln.- fl 2 - — -- + ^ + - + — 

and so on. 

The control function to(g) is defined from the quasifixed-point condition 



|W) = 0, 



which gives 

u(9) = 



Q \ V2 



This control function in the weak-coupling limit, as g — > 0, behaves as 



and in the strong coupling limit, as g — > 00, it has the asymptotic behaviour 

For g e [0, 00), function (69) changes in the interval \J~2 < oj{g) < 00. 
Defining 

fk(g) = F k {g,w{g)), 
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from (66) and (67) we have 



Similarly, calculating F2(g,ui), we come in the second order to 

/2(9)= W-i-i-^ + ii. (74) 

Wishing to estimate the accuracy of the approximations fi(g) and /2(g), let us notice that the exact function (61) 
changes from —00 as g — ► to +00 as g — > 00, crossing zero at g = g c = 2.758, that is, f(g c ) — 0. Therefore, we 
cannot define the percentage error in the standard way as 100% x [fk(g) — f{g)]/f{g), since such a definition contains 
zero in the denominator. Instead of this, we may evaluate the accuracy of a crossing point given by the corresponding 

(k) 

approximation, that is, the accuracy of the solution g c to the equation 

f k (gi k) ) = 0. (75) 
For the first approximation, we have — 0.585 which gives the error of —79%, as compared to the exact g c — 2.758. 

(2) 

For the second approximation, we get gc — 1.352, whose error is —51%. As is seen, this accuracy is not high, so 
that it is desirable to improve it. 
Introduce the function 



i yr+%-1 

a{9) = = 3. 9 • (76) 



This function changes in the interval < a(g) < |, with the asymptotic behaviour 

1 3 

a ^~2^8 9 fa - * )' 

a(g) ~ — ^= - J- (g — » oo). 
V¥ 3. 9 vy y 

Using (76), we can write for the approximations fk{g) the following expressions: in the zero order, 

*<»)— 5 h "-s^y (77> 



where a = a{g), in the first order, 



and in the second order, 



h(g) = fo(g)-\ + \a, (78) 



/2(5H/o(s)-^(l + 5a-14a 2 ). (79) 

The zero-order approximation (77) correctly describes the weak and strong coupling limits of the exact function 
(61), but it is not accurate in the intermediate region. In this region, the behaviour of the approximation 

fk(g) = fo(g)+Pk(a) 

is governed by a series Pk{a) in powers of a = a(g). Renormalizing this series twice, according to the bootstrap 
procedure [16], we obtain 



fZ{g) = fo(g) - ^ cx P |5acx P } 



(80) 
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where a = a(g) is given by (76). The obtained self-similar approximant (80) provides much better approximation to 
function (61), compared to fi(g) and /2(g)- The crossing point g* = 2.858, defined by the condition f^igl) = 0, is 
quite close to the exact g c and gives an error 3.6%. That Eq.(80) is an accurate approximant is also clearly seen in 
Fig. 1. 

B. One Dimensional Oscillator 

From the model of the previous section, we now pass to a more realistic case of the double-well oscillator with the 
Hamiltonian 

" = 4;£ + i-r 2 + ^ < 81 » 

in which x G (—00, +00) and g G [0, 00). The problem of finding the eigenvalues of Hamiltonian (81) is a challenge for 
any analytical method, although there are several numerical techniques calculating the eigenvalues with reasonable 
accuracy [38-43]. It is especially difficult to calculate the lowest energy levels. The main problem here is that instanton 
contributions are crucial in the weak-coupling region providing for an exponentially small splitting of energy levels. 
In addition, the energy of the ground-state level is not a monotonous function of the coupling parameter g, which is 
also related to the instanton contributions. Below we shall consider the most difficult case of the ground-state energy 
and that of the first excited level separated from the former, in the weak-coupling region, by an exponentially small 
gap. 

To construct interpolation formulas, we need asymptotic expansions for the weak and strong coupling limits. We 
shall use such expansions derived in Ref.[44]. The ground-state energy e+{g) corresponds to a symmetric wave 
function, while the first excited level, with an energy e_ (gr) , corresponds to an antisymmetric wave function. These 
energies can be written in the form 

e±(g) = e(g) T ^(g), (82) 

in which 

e{g) = \[e + {g)+e_{g)], A(g) = e_(g) - e+{g). (83) 

Let us notice that the Hamiltonian (81) is shifted, as compared to the standard form, by the term 1/16(/, which makes 
the spectrum of (81) everywhere positive [44]. 

For the average energy and the gap, defined in (83), we have [44] in the weak-coupling limit, when g — > 0, 

1 21 

~e(9) ^ T2 ~ Va 9 (84) 

and, respectively, 

AfoJ-^expQ + c), (85) 

where 

303 , V2 9 

a = , b = , c = — . 

1024' 4 ' 4 

In the strong-coupling limit, when g — > 00, we may derive [44] for the energies 

e + (g) * |(6 fl )V3 - I (65) -V3 + l| (6fl) -i _ §g(6 5 )- 5 /3, 



e-(g) - |(10)V" - |(10 ff )-V» + g!(10 5) - - ^dOg)-^. 
From here, for the average energy we get 

e(g) ~ V /3 + Bg- 1 ^ + CcT 1 + DcT 5/3 , (86) 
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where 



A = 1.552580, B = —0.242850, C = 0.221181, D = -0.074868. 
And for the gap, we find 

A(g) ~ A 1 g 1 / 3 + B^ 1 ' 3 + C^ 1 + D.g- 5 / 3 , 

where 

Ai = 1.742319, Bi = -0.210539, d = 0.081250, D 1 = -0.070721. 
Constructing a self-similar approximation from the right to the left, we have the right approximant 

e* 4 (g,oo) = e*( 3 )T^A*( 3 ), 



(87) 



in which 



e*(9) = 



i + A 4 g 4/3 exp 



AB 
Ag 2 / 3 



-,1/4 



Dg 



g W + (%\D\) 



1/4' 



(88) 



(89) 



The most difficult here is to interpolate between the power-law expansion (87) for the gap, in the strong coupling 
limit, and the exponential behaviour (85) in the weak-coupling region. Nevertheless, employing the technique of Sec. 
II, we obtain for the gap the form 



A* (ff ) = ^W^ 



(90) 



describing a renormalized instanton contribution, where 



«* (.9) = 



3/2 , A 3/2 2 / 3-Bl 

a ' + A{ f exp < — — - 



\2A l9 3 ' 2 6XP \B l9 V 3 



2/3 



0*{g) = T = 9 h^T - 0.22417. 

[|£i| 4 /3 + (|6|^i) 4 / 3 (r + g 2 / 3 y/ 2 gf \ A i\ b \ J 

The behaviour of two branches of Eq. (88), compared with the exact numerical data, correctly describes the 
nonmonotonic behaviour of the ground-state energy and the exponential branching at g = 0. The accuracy of 
the constructed formulas is very good for both the instanton-dominated region (g <C 0.2) and instanton-free region 
(g 3> 0.2), with the error tending to zero in both these limits. The most difficult for description is a narrow intermediate 
region around the point g w 0.2, where an error is about 25%. The accuracy can be improved taking into account 
more terms of asymptotic expansions, but then the formulas become essentially more complicated. 

C. Quasistationary States 

Quasistationary or resonance states are encountered in a variety of studies in atomic and molecular physics. A 
good discussion and many references are given in Rcf. [45]. Several numerical calculations are known for this 
problem [6,42,45,46]. Here we derive analytical formulas for both the real and imaginary parts of the spectrum of the 
Hamiltonian 

"—tst + t*-*- < 91 » 

in which x G (— oo, +oo) and g e [0, oo). 

First, we have to analyse the asymptotic behaviour of the spectrum in the weak and strong coupling limits. To this 
end, we invoke perturbation theory starting from the Hamiltonian 

1 d 2 u 2 9 , . 

Fo = -2^ + T*< ^ 
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in which u is a trial parameter to be converted into a control function. Introducing, for convenience, the notation 

E k (g,u)= (n+^j F k (g,u) (93) 

of a /c-order perturbative expression for the energy of a level n = 0, 1, 2, . . ., we find 

u 

F (g,u) = u, Fx{g,u) =u- -(2a- (3), 

F 2 (g,u) = F^u) -^(a 2 - 2al3 + 2/5 2 <5), (94) 

and so on, analogously to the case of the anharmonic oscillator [13,31], with the notation 

1 a _ 675 _n 2 +n+ 1/2 _ 17n 2 + 17 n + 21 

a=l 2 ' P = 3"' T= rTTo — ' " = • 

n+1/2 (67) 

The control function u = u(#) can be defined from the fixed-point condition 

^-F k (g,u) = Q. (95) 



In the first order, we get the equation 



as a result of which 



u A - u + 675 = 0, (96) 



a = f3= — . 

u z 

Substituting the solution u(g) to Eq. (96) into (94), we define 

h(g) = F k (g,u(g)). (97) 

Then, from (94) we have 

h{9) = \ u + h h{9) = h{9) + ^ (1 ~ 25)a2u - (98) 

From three solutions of Eq. (96) we need to choose that which satisfies the boundary condition u(g) — ► 1, as g — * 0. 
Such a solution is 

u(g) = Q v / 729 7 2 .9 2 - 3 - 3 7ff ) ' + l - Q v /729 7 2 3 2 - 3 - 3 7 fl) ' ■ (99) 
The control function (99) is real for g < g n , where 

9n = 30 = n + 1/2 v| = 

y 7 n 2 +n+l/2 y ' y 27 v 7 

Complex roots of (99) appear only after g > g n . The fact that Im u k (g) = for g < g n leads to Im f k (g) = 0, when 
9 < g n - 

Asymptotic expansions in the weak and strong coupling limits can be written for arbitrary energy levels. For 
illustrative purpose, we shall write down expansions for the ground state (n = 0) and the first excited state (n = 1). 
For the ground state, when n = and 7 = 1, function (98) yields for the real parts 

^ h(g)^i-\-\g 2 - vg* - ™g\ 

Re Mg) *\-\g- ^ - ^g 3 - ^ (n = 0), (101) 
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while the imaginary parts are 

Im h(g) = Im f 2 (g) = 0, g -> 0. (102) 

For the first excited level, for which n = 1 and 7 = 5/3, for the real parts of (98) we find 

, , s 5 25 o o 13125 4 

Rc ACg) ~ 1 - - —g 2 - 125.g 3 - ~^g\ 

Re/ 2 ( 3 )~l-- 3 --g 2 -— g 3 -— g 4 (n = 1), (103) 

while their imaginary parts are again as in (102). 

In the strong-coupling limit, when g — > 00, for the ground state we obtain the real parts 

Re h {g) * l(6g)^ + \(&g)-V* + y^g)^ + y^)" 573 - ^r V \ 

Rc f 2 fg) ~ -(Gg) 1 / 3 + ^(Gs)- 1 / 3 + —(eg)" 1 + ^-(6g)- 5 / 3 - J^l( 6g )-7/3, (104) 
J w; 96 v y; 288 v y ' 144 v y> 1944 v y; 23328 v ' v; v y 

and the imaginary parts 

to m a -2^! (%) V3 + ^ (%) -V3 _ ^ (%) -V3 _ ^ (69 ,-V3. (105) 
For the strong-coupling limit, in the case of the first excited level (n = 1), we find the real parts 
Re Mg) ~ ^(Wg) 1 / 3 + i(lO.g)- 1 / 3 + ^(lO^)" 1 + ^(lO*?)" 573 - ^(lO*)" 7 ", 
Re /„(*) * ^(lO.g) 1 / 3 + |(10.g)-/3 + | (10fl) -i + ^(i 0ff )-«/3 _ ^(lO,)"^, (106) 
and imaginary parts 

Im / l( .9) * -^(10 fl )V3 + ^(I0.g)-V3 _ ^(lO,)" 5 / 3 - (10. 9 )- 7 /3, 

Im / 2 ( 3 ) - -^(10. 9 ) 1/3 + ^(10^)- 1/3 - W (105) " 5/3 - ^l (105) " 7/3 - (1 ° 7) 

The accuracy of the real parts of the approximations in (98) is sufficiently good; the maximal error in the first order 
is —3.3% and that of the second order is 2%. In the case of the corresponding imaginary parts, the maximal error 
for g > g n is on the order of 10%. However, for g < g n , imaginary parts are identically zero, because of which their 
weak-coupling expansions do not exist. 

To correct the described deficiency, let us consider the fixed-point condition (95) of second order. Then, the 
equation for the control function becomes of the form 

u 6 - 2u 4 + 16 7 .gu 3 + u 2 - 16751* + 120j 2 Sg 2 = 0. (108) 

For the ground-state level, for which n = 0, 7=1, and S = 7/12, Eq. (108) reduces to 

u e - 2u 4 + I6gu 3 + u 2 - IQgu + 70.g 2 = 0. (109) 

The solution to (109) will be denoted by u*(g), in order to distinguish it from that of (96). Among the solutions 
of Eq. (109) it is necessary to choose that which satisfies the asymptotic boundary condition u*(g) — ► 1, as g — ► 0. 
Substituting u*(g) into 

p, x 3 3 3g 1 3g 21g 2 
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we come to 

f^g) = F 2 (g,u*(g)). (110) 

Defining the real and imaginary parts of (110), we take those of them for which 

Re/ fc *( 5 )>0, ImA*( S )<0. 

To solve Eq.(109), we introduce a new variable 

A = u(u 2 -1), (111) 

for which (109) reduces to a much simpler equation 

A 2 + I6g\ + 70<? 2 = 0. (112) 

From two solutions to this equation, Ai ; 2 = — (8 ± iV&)g, we need to choose such that, together with Eq. (Ill), gives 
the control function u*(g) satisfying the conditions discussed above. The desired solution to (112) is 

A(s) = -(8 + iv%. (113) 
Then from the equation u 3 — u — A = we find 

u*(g) = Q + ^VS^) + | Q + , (114) 

where A = X(g) is defined in (113). 

In the weak-coupling limit, when g — > together with A — > 0, approximant (110) has the following asymptotic 
expansions for the real part 

Re «W * 1 - \ S - 2 i<S - 41/ - ifV - 5547/ - ^ 9 «, (in) 
and for the imaginary part 

Im /!(,) * -\^g* - 27V6. 9 4 - ^V6 5 5 - ^V6 5 6 . (116) 
In the strong-coupling limit, when g — > oo and |A| — > oo, for the real part of (110), wc have 
Re fi(g) ~ 0.672436 3 1/3 + 0.145202. 9 ~ 1/3 + 0.016735.9" 1 + 0.000603<T 5/3 - 

- 0.000088cT 7/3 - O.OOOOlO.g" 3 - 0.7 x lO^V 11 / 3 , (117) 
and for its imaginary part, we find 

Im fi(g) ~ -1.155930 5 1/3 + 0.246594<T 1/3 - 0.000750,T 1 - 0.001364 ff - 5/3 - 

- 0.000101.g- 7/3 + 0.000003 ff ~ 3 + 0.000001<T 11/3 - (118) 

The approximant (110) possesses the weak-coupling expansion for its imaginary part, thus, correcting the deficiency 
of the approximations in (98). The values of the real and imaginary parts of f£(g) for different g, as compared to the 
precise numerical calculations [42,47], are given in Table III. The maximal percentage error of the real part is 0.7%. 
The percentage error of the imaginary part cannot be correctly defined, since Im f(g) — > as g — > 0. Recall that the 
ground-state energy E(g) for Hamiltonian (91) is related, according to (93), with the function f(g) as 

E{g) = \f{g). (119) 

To write an analytical formula approximating f(g) in the whole range of the parameter g G [0,oo), we may use 
the asymptotic expansions derived above. Employing the technique of Sec. II, we can obtain for the real part the 
approximants 
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Re / 2 *(o, 0) = [exp(-6oi</) + A\g 



6 21 1/6 



(120) 



Re / 3 *(o,0) = [cxp{-6( ai + a 2 ) 3 } + A 6 l9 2 ] 1/6 + a 2 g 



. 3/2 



2/3 



(121) 



in which 



ai = 1.5, a 2 = 0.548, A x = 0.672436, A 2 = 0.145202. 



The behaviour of these approximants is shown in Fig. 2. The maximal percentage error of f 2 {g, 0) is 24% and that 
of /!(<?, 0) is 9%. 

To reach good accuracy for the imaginary part, we need to consider higher-order approximations. For the imaginary 
part we find the following sequence of self-similar approximants: 



exp 



36 \ fa 
Ta 9 ) + U 



ImfZ(g,0) = -ag 3 
where a = |\/6, 6 = 27 VE, and A = 1.155930, 

lmf* i (g,0) = -ag 3 
where a, b, and A are the same as in (122) and B\ = 0.477045 

816 \ 81 



") 3/ V 



-4/3 



276 \ 27 n 2 



7/6 



9/8 



-8/9 



lm/*(g,0) = -ag 6 



CXP l"7b^ 5 J + 70^ 2 



7/6 



27 3 
20^ 



-i 10/9 



3/2 



-2/3 



where B 2 = 0.178716 and B 3 = 0.496502, 
Im /*(<?, 0) = -ao 3 



exp 



2436 \ 243 



182a 182 



7/6 



52 5y 



-1 10/9 



+ 



"26^ 



13/12 



+ 



(3) 



15/8 



-8/15 



(122) 



(123) 



(124) 



(125) 



where B 4 = 0.073553, B 5 = 0.196401, and B 6 = 0.552922. 

The behaviour of the self-similar approximants (122)-(125) is displayed in Fig. 3. The curves corresponding to 
Eqs. (124) and (125) almost coincide in this picture. Fig. 3 clearly demonstrates the convergence of the sequence 
{Im f£(g, 0)} to exact data marked by diamonds. 



V. NONLINEAR HAMILTONIANS 

Here we show that our approach is applicable not only to linear problems of quantum mechanics but to nonlinear 
problems as well. 

A. One Dimensional Case 

Consider the nonlinear Hamiltonian 



H 



1 d 2 1 



x 2 +gip 2 (x), 



2dx 2 2" 

in which x € (—00, +00), g G (—00, +00), and ip(x) is a wave function. 



(126) 
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The Hamiltonian (126) is a prototype of the if 4 model of quantum field theory. There exists a controversy in 
the interpretation of the so-called "triviality" of this theory (see discussion in Refs. [48-50]). Therefore, developing 
methods that could successfully deal with nonlinear Hamiltonians of the type (126) could be useful for (p 4 quantum 
field theories, as well as for those field theories that include scalar-field terms, like (126), in their Lagrangians, e.g., as 
in the Higgs model [51]. Another important application of the nonlinear Hamiltonian of the form (126) is for describing 
properties of atoms confined in magnetic traps [52-54]. Such magnetically trapped atoms of 87 Rb, 23 Na, and 7 Li, 
as has been observed recently [55-57], can exhibit the phenomenon of Bosc Einstein condensation. The possibility of 
the direct observation of Bose condensation distinguishes these alkali gases from liquid 4 He where this condensation 
could be investigated only indirectly (see the related discussion in Refs. [58-63]). The system of condensed trapped 
atoms corresponds to the ground state of the Hamiltonian (126). We consider here the one-dimensional case which 
serves as an illustration of the applicability of the method. 

It is worth emphasizing that the coupling parameter g in Eq. (126), when a system of TV condensed atoms 
is considered, is proportional to N. Because of this, for N 3> 1, the coupling parameter can become large and, 
consequently, perturbation theory in powers of g is of no sense. Although some thermodynamic characteristics of 
trapped atoms can be approximately analyzed disregarding their interactions, that is, nonlinearity [64-67], but a 
correct description certainly has to take into account these interactions [68], i.e., nonlinearity, since for the cases 
related to experiment [55-57] the effective coupling is strong, j > 1. 

To derive analytical formulas for the spectrum of the Hamiltonian (126), we need to find appropriate asymptotic 
expansions. Let us start with the Hamiltonian 

1 d 2 v 2 

H « = 4^ + T* 2 < < 127 > 

in which u is a trial parameter. Employing perturbation theory with respect to the perturbation 

1 



we may find the fc-order approximation 



AH=^(l-u 2 ) x 2 +g^ 2 (x), 



E k (g,u)= (n+^jF k {g,u) (128) 
for the energy levels labelled by n — 0, 1, 2, . . .. For the function F k defined in (128) we have 

F (g,u)=u, F 1 (g,u) = ^u+^j+—^gV^, (129) 



where 



1 f +cc 



H n (x) being a Hermite polynomial. In particular, 



\/2V 4\/2V 64^ 

The control function u = u(g) can be found from the fixed-point condition 

d 



. )u F 1 ( fl) u) = 0. (130) 



The latter, with the notation 



a=—^- gi (131) 

n+l/2 y ' v ; 

yields 

u 2 + u 3/2 a - 1 = 0. (132) 
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Substituting the solution u(g) to Eq. (132) into (129), we get 

fk(g) = F k {g,u{g)). 

For instance, 

h(9) = \ ~uj , u = u{g). 
In the weak-coupling limit, a — > when g — > 0, according to (131). Then (134) gives 



(133) 



(134) 



fi(g)^l + a 



\ a + 32° -US a + 2048 



a 5 . 



In the strong-coupling limit, when g — > oo together with a — > oo, Eq. (134) leads to 

/ 1 ( 3 )4« 2/3 + ^- 2/3 -^- 2 + ^- io/3 . 



(135) 



(136) 



In the case of a negative coupling parameter g < 0, the weak-coupling limit, when 5 — > — and a — ► —0, gives the 
same asymptotic expansion as (135). However, the strong-coupling limit, when 5 — > —00 and a — > —00, is different 
from (136). If a — » —00, then Eq. (134) behaves as 



cr 1 1 3 13 

h{9) - - y + 2a*~2^ + 2a™ ~ 2a^' 



(137) 



In the region of negative g < 0, the energy (128) is positive for small \g\, and, as g diminishes, the energy becomes 
zero at a critical value g c . The latter can be found from the definition 



A(.9c) = 0. 

The form (134) shows that equality (138) holds true for v? c = 3. Then, Eq. (132) immediately gives 

2 



Because of the relation (131), one has 



OL c = - 



9c 



3 3 / 4 



-0.87738. 



2 Jn 



For example, for the ground-state level, with n = 0, one finds 

& = -1.09964 (n = 0). 



(138) 



(139) 



(140) 



(141) 



For g < g c , the energy becomes negative, which implies the instability of the system. 

Note that invoking the notation (131), we have managed to write the asymptotic expansions (135), (136), and 
(137) for arbitrary energy levels. Using the derived expansions, we can construct analytical formulas for the whole 
spectrum of the considered nonlinear Hamiltonian. Thus, following the standard way of Sec. II, we find for a > the 
approximant 



/ 2 >,0) 



1 + axAlJ&Ala 



2/3 



Afa 2 



1/3 



(142) 



where a\ = 1, A\ = 3/2, and A 2 = 1/2. Similarly, we can write down 



/ 2 * (a, 00) = 

where A\ and A 2 are the same as before and 



5/4' 



1/4 



(143) 



T 2 = 



A 2 \2 12 A{ 



at 

12 411 



1/5 



= 0.776. 
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For a < 0, we may construct 

/|(a, -0) = 1 + a x a 
with a 2 = —1/8 and B\ = —1/2. Another approximant is 

/ 3 >,-oo) =B ia 2 + S 2 
where B 2 = 1/2, B 3 = -1/2, and 



/ 2a 2 \ 


B-. 9 


cxp a 


+ -k-OL 


\ a i J 


a\ J 



1/2 



Bi + \a\ [a 



16B 3 



-r 3 



7/4' 



-1/4 



(144) 



(145) 



7 /af B 9 2 3 xi/7 



73 



I-B3 



t l x -'2 
2 20 



0.198. 



The behaviour of approximants (142)-(145) is such that /2*(a,0) practically coincides with /|(a,oo) for a > 0, 
while /I (a, — 0) coincides with /|(a, — 00) for a < 0. A slight difference between these approximants is noticeable 
only in the region — 1 < a < 1, as is shown in Fig. 4. 

B. Radial Model 

In the previous section we considered a one-dimensional nonlinear problem. It is straightforward to apply the same 
approach to nonlinear problems of higher dimensionalities. As an illustration, we consider a spherically symmetric 
model with the nonlinear radial Hamiltonian 



1 d 2 1(1 + 1) 1 2 n2 / x 



(146) 



in which r G [0, 00), n = 0, 1, 2, . . . , I = 0, 1, 2, . . ., the coupling parameter g e (—00, +00), and R n i is a radial wave 
function. 

Starting from the trial Hamiltonian 



1 m + i) v? 2 
^ = -2^ + ^^ + T r ' 



(147) 



with a trial parameter w, we invoke perturbation theory with respect to the perturbation AH = H — H - For the 
fc-order approximation of the spectrum we may write 



E k {g,u)= \ 2n + l + -)F k {g,u). 



Introduce the effective coupling parameter 



a = 



Jnl 



2n + / + 3/2 



9, 



(148) 



(149) 



in which 



>nl 



2nl 



r(n + / + 3/2) 



r 4(;+i) e -2r^ 



L l n +1/2 (r 



1 4 



dr, 



where T is a gamma-function and L l n is an associated Laguerre polynomial. For the particular cases that will be 
analysed in what follows, we have 



^00 = 



-An = 



35 



J10 = 



147 



J02 = 



231 



2V2tt 24V2tt 128V27T 160\/2tt 

Accomplishing the same steps as in the previous section, we can find the spectrum 



e nl (g) =[2n + l + -) A( 3 ), 



(150) 
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in which fi(g) is defined identically to (133). The control function u(g) is given by the solution of equation (132), 
only with a introduced in (149). With this renotation for the parameter a, all asymptotic expansions for fi(g) are 
the same as in Sec. VII. Therefore, the corresponding approximants /£ will have the same forms (142)-(145). 

What differs the considered nonlinear radial model from the one-dimensional nonlinear case is that the energy 
levels are labelled now by two quantum indices, n and I. For some combinations of these indices, the specific effect of 
level crossing, when varying g, may occur. To illustrate this, let us analyse several lower levels of the spectrum (150). 
In the weak-coupling limit we may find 

3/ g \ 5/ 7g \ 

eoo "2( 1+ v^J' 601 " 2 v + ) ' 



e 10 ~-(l + 21 t- ) , e 2^-(l+ 33 !_ ) , g^O. 

While for the strong-coupling limit, we have 

eoo(g) 1.219334 5 2 / 3 , e w (g) ~ 1.353502 5 2 / 3 



eoi(g) 1.418783 5 2 / 3 , e 02 (g) ~ 1.576577 5 2 / 3 , .g -» oo. 

As is seen, at small g the energy level eoi is lower than eio, but at large g, vice versa, the level eio becomes lower 
than eoi- 

Another difference with the one-dimensional nonlinear case is the value of the critical coupling constant g c at which 
the energy level e n i(g) crosses zero. This critical parameter is defined as in (138) and gives the same a c as in (139). 
But, because of the different relation (149) between a and g, we now get, instead of (140), the critical coupling 

2n + l + l)X l 

For the ground state-level, instead of (141), we find 

g c = -2.19927 (n = I = 0), 
so that | gc | is about twice larger than that for the one-dimensional case. 



VI. WAVE FUNCTIONS 



In this section we show that the approach we have developed permits us to construct analytical expressions not 
only for energy levels but for wave functions as well. We shall concentrate our attention on the most interesting, from 
our point of view, case of nonlinear equations. 

A. Nonlinear Schrodinger Equation 

Consider the equation 

~l^ + l x2 ^ + gi,3 = E ^ (151) 

in which x G (— oo, +oo) and g G [0, oo). The nonlinear Schrodinger equations of this type are used for describing 
Bose-condensed atoms in magnetic traps [55-57]. This kind of equations is also often called the Gross-Ginzburg- 
Pitaevskii equation [69-72]. Condensed atoms correspond to the ground state of (151), which is assumed in what 
follows. 

The wave function ip = ip(x) is normalized by the condition 

/+oo 
\i/>(x)\ 2 dx = l. (152) 
-oo 

From Eq. (151) it follows that the wave function is an even function, so that ip(x) — tjj(-x). Because of this, it has 
an expansion 
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tp(x) ~ c + c 2 x 2 + c 4 x 4 , (153) 

as x — * 0, in even powers of x. For large x, the harmonic term in (151) becomes predominant, and the wave function 
has the asymptotic behaviour 

ip(x) ~ Aexp (-^ 2 ) ( 154 ) 

as x — > 00. Our aim is to find an analytical expression for ip{x) valid in the whole region of x g (—00, +00). Note 
that Pade approximants are not able to interpolate between such different types of behaviour as the power law in 
(153) and exponential in (154). But using the self-similar interpolation of Sec. II, we easily obtain the self-similar 
approximant 

tp*(x) = Acxp j-ix 2 + ax 2 cxp(-fex 2 )| , (155) 

where ip*(x) = ip^(x,0) and 



A = c , a=H , b= 2- 

2 c acg 



x 1 



Expression (155) acquires a transparent physical meaning when written in the form 

ip*(x) — Acxp - 

in which 



£,{x) = Q - acxp{-6x 2 }^ 



-1/2 



plays the role of an effective correlation length. 

The parameters a, b, and A in (155) are not independent. The relation between them can be found if we expand 
(155) in powers of x and substitute this expansion into Eq. (151). For small x, Eq. (155) gives 

i>*(x) — Co + c 2 x 2 + C4X 4 (x — > 0), 

where 64 = C4 + c 2 /2co- Substituting this into (151) and equating the terms at like powers of x, we find the relations 

c 2 = c (gc 2 -E), c 4 = ^ (1 + AgclE - AE 2 ) . (156) 



The latter yield the equalities 



1 ~ 2(l-2aVE-l 
a=- 2+9 A 2 -E, b = J A (157) 



showing that among four parameters, a, b, A, and E, there are only two independent. Two additional equations for 
defining all parameters are the normalization condition (152) and the definition of the energy 

E*(g) = (r,Hr), (158) 

where the Hamiltonian H is the same as in (126) and the notation E*(g) = E stresses that the energy is obtained by 
using the self-similar approximant (155). 

The values of a, b and E* depend on the coupling parameter g. Thus, for the weak-coupling limit j^Owe have 



Then, relations (156) give 



c 2 _ 1 y/2-1 c 4 _ y/2-1 
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Respectively, from (157) we get 



a~ ^ 1 g = 0.165247ff, 6 ~ = 0.235702. 
V27T 3y2 



This demonstrates that the function (155) reduces to the Gaussian form when g — > 0. 
The variational Gaussian function 

^)=Q 1/4 cxp(-Jr 2 ) (159) 

is often used not only for small g <C 1 but for arbitrary g £ [0, oo), with the effective frequency u = u(g) defined by 
the minimum of the energy {tpc Hipa), which gives 

u 2 +u 3 / 2 ^-g-l = 0. 

Such a variational energy is very close, with a deviation not more than several percent, to the energy 

Z V 7T 



corresponding to Eq. (142), which results in 



27 \ 2/3 . 27 2 



9 + 1T9 



4v/2tt V 4tt 



1/3 

(160) 



Another very often used approximation for treating Bose-condensed atoms in harmonic traps is the Thomas-Fermi 
approximation, (see e.g. [73-75]) which for Eq. (151) leads to the wave function 

/ x 2 _ x 2\ V2 
■>Ptf(x) = I — — — J , \x\ < x , 

iPtf(x) = 0, \x\ > x , (161) 

in which 

3 ^ 1/3 



x a = \JlE T F = {^9 

The energy in the Thomas-Fermi approximation is obtained from the normalization condition (152) giving 

1 /3 \ 2/3 

E TF (g) = -[-g) . (162) 

The Thomas-Fermi approximation is assumed to be valid for large g — > oo. However, even then the wave function 
(161) is correct only for x <C x , where it has an expansion 

iPtf(x) ~ c + c 2 x 2 + c 4 a; 4 

with the coefficients 

x c c 

The behaviour of the function (161) near the boundary x = x is not correct. Also, this function is not appropriate to 
evaluate the mean kinetic energy, producing a divergence for any g (see discussion in [73,75]). In order to understand 
when the total energy (162) for the Thomas-Fermi approximation starts giving reasonable results, we present in Fig. 
5 the energies (158), (160), and (162). The first two energies, E*{g) and E^g), almost coincide with each other, 
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having correct asymptotic behaviour in the weak as well as in the strong coupling limits. The Thomas-Fermi energy 
ETF(g) possesses an incorrect weak-coupling limit and becomes a reasonable approximation starting from g = 7. 
The density 

n (x) = \ip{x)\ 2 (163) 

for the corresponding wave functions and g = 0.2 is presented in Fig. 6, where the density n* (x) — | ip* (x) | 2 of the self- 
similar approximation (155) practically coincides with the density no(x) = \ipa{x)\ 2 of the Gaussian approximation 
(159), as it should be in the weak-coupling limit. In this limit, the behaviour of the density riTF(x) — \iPtf(x)\ 2 of 
the Thomas-Fermi approximation is not correct. As is known [73,75], the latter approximation is incorrect near the 
boundary even in the strong-coupling case. Then the density n(x) in the self-similar approximation is close, except 
near the boundary, to that of the Thomas-Fermi approximation. The self-similar approximation n*(x) coincides with 
nTF{x) for small x and smoothes the incorrect behaviour of utf(x) around the boundary. In the strong-coupling 
limit, the density nc(x) of the Gaussian approximation is not accurate. 

The direct evaluation of the accuracy of each approximation can be done by calculating the residual term 

R(x) = Hip(x) - (ip, Hip) (164) 

for Eq. (151), where H is defined in (126) and ip(x) is a wave function of the corresponding approximation. The 
residual terms for g 3> 1 for the self-similar approximation (155) is practically zero, meaning that (155) is an almost 
exact solution of Eq. (151). For the Gaussian approximation (159), the residual term is much larger, telling that 
this approximation is much less accurate. And the residual for the Thomas-Fermi approximation is divergent at the 
boundary point xq, though far from this point it is close to zero. 

The integral characteristic of accuracy of the corresponding solutions is the dispersion 



a(V0 



\R(x)\ 2 dx 



1/2 



(165) 



We calculated this quantity for < g < 100. The maximal dispersion for the self-similar approximation (155) is 
around 1, for the Gaussian approximation (159), it is about 20, and for the Thomas-Fermi approximation, it is 
divergent. 

In this way, the self-similar wave function (155) is the most accurate solution of the nonlinear Schrodinger equation 
(151), as compared to the Gaussian and Thomas-Fermi approximations. This function (155) represents the exact 
solution very well for all x and g. In the weak-coupling limit g — > 0, it becomes close to the Gaussian form, and 
in the strong-coupling limit, it approaches the Thomas-Fermi wave function for all x except the boundary where it 
smoothes the incorrect behaviour of the latter function. The crossover point between the weak-coupling and strong 
coupling regimes occurs, as numerical calculations show, at around g s=s 5. This crossover point can also be evaluated, 
by an order of magn itude, analytically as follows. Notice that the characteristic length for the Gaussian function 
(159) is xq — y/2/u with u w 1, and that such a length for the Thomas-Fermi function (161) is xq — (Sg/2) 1 ^ 3 . 
These characteristic lengths, typical of the weak-coupling and strong-coupling regimes, respectively, coincide, that is 
x G = x Q , at g w 2 5 / 2 /3 w 2. 

B. Vortex Filament Equation 

Now we shall show that our approach permits us to find accurate analytical approximations for the function 
describing the structure of vortex filaments. Considering an unbounded Bose system and making in the nonlinear 
Schrodinger equation the substitution ip(r) = f(r)e t(f , in which r and ip are dimensionless polar coordinates, one 
comes [70,72] to the equation 

g + ;f-i + /-/ 3 = »- (i«> 

The solution to this equation is usually obtained numerically [70,76,77]. Here we shall construct a sequence of 
analytical approximations for the solution to Eq. (166) and compare them with the known numerical data. Note 
that the equations similar to (166) have been considered as well for describing magnetic solitons [78], isomeric states 
of quantum fields [79], and vortices of complex scalar fields [80]. Therefore, the possibility of deriving accurate 
analytical solutions to these equations is important for many applications, such as condensed Bose gas, superfluid 
helium, magnets in strong magnetic fields, and different models of quantum fields. 
At small r — > 0, the solution to Eq. (166) has the asymptotic expansion 
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f(r) ~cr(l- ^ , (167) 

where c is a constant. At large r — > oo, one gets 

f(r) - 1 - \r~ 2 - jjr-* - ^r" 6 . (168) 
Employing the approach of Sec. II, we easily obtain the following sequence of self-similar approximants: 



J'l '('■■»)= r 2 r(l + ir 2 



-1/2 



-1/4 



/ 4 *(r,0)= C4 r(l + i r 2 + -r 4 + -r 6 



-1/6 



-1/8 



/b (r, 0) - c 5 r + ^ + + (169) 

in which the coefficients, defined so that to give the correct asymptotic expansions, are 

c 2 = 4~ 1/2 = 0.5, c 3 = 4~ 1/4 = 0.707, c 4 = 16~ 1/6 = 0.630, c 5 = 140~ 1/8 = 0.539. (170) 

The behaviour of the approximants f%(r, 0) is shown in Fig. 7, compared with numerical data [70,76,77]. As it can 
be concluded from this figure, f^(r, 0) is a very accurate solution. 

VII. CONCLUSION 

We have developed an approach for obtaining analytical solutions of quantum-mechanical problems. This approach 
makes it possible, starting from asymptotic expansions having sense only in the vicinity of limiting points, to derive 
interpolation formulas valid in the whole range of variables. The developed method is rather general and can be ap- 
plied to various problems. We demonstrated its applicability to several quantum mechanical models, such as different 
anharmonic oscillators, double-well potentials, resonance models with quasistationary states, and nonlinear Hamil- 
tonians. The method permits one to construct accurate analytical expressions for energy levels as well as for wave 
functions. It is important that this method provides a regular procedure for deriving a convergent sequence of subse- 
quent approximations, so that it is possible to reach the desired accuracy by calculating higher-order approximations. 
The idea of the approach is based on the self-similar approximation theory [8-17], this is why we call the method 
developed in the present paper the self-similar interpolation. The method can find numerous practical applications, 
for example, for analysing spectral properties of atoms and molecules, for studying the physics of quantum dots, and 
for investigating the behaviour of Bosc condensed gases in magnetic traps. 
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TABLE I. Percentage errors of the self-similar approximants for the ground-state energy of the sextic oscillator, as compared 
to numerical data for e(g). 



3 


e(9) 


£i (3,0) 


£i(3,°o) 




£2(3,0) 


£3(3,0) 


£4(3,0) 


£5(3,0) 


(3,0) 


0.1 


0.586945 


-8.25 


7.25 


-0.50 


-5.45 


-3.99 


-3.06 


- 2.46 


-1.95 


0.5 


0.717813 


-10.56 


5.88 


-2.34 


-4.71 


-2.48 


-1.44 


- 0.88 


-0.53 


2 


0.915219 


-8.47 


3.99 


-2.24 


-2.40 


-0.82 


-0.33 


- 0.11 


-0.06 


50 


1.858487 


-2.47 


1.03 


-0.67 


-0.16 


-0.03 











1000 


3.850896 


-0.59 


0.25 


-0.16 


-0.01 
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TABLE II. Percentage errors of the self-similar approximants for the ground-state energy of the octic oscillator. 
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e(9) 


e! (9,0) 


£1(9,00) 


4(g) 


£2(9,0) 


£3(9,0) 


£4(9,0) 


4(9,0) 


4(9,0) 


0.1 


0.620514 


-12.01 


7.82 


-2.10 


-7.90 


-5.72 


-4.27 


-3.30 


-2.66 


0.5 


0.745510 


-12.54 


6.10 


-3.22 


-5.77 


-3.15 


-1.88 


-1.21 


-0.81 


2 


0.911090 


-9.66 


4.39 


-2.64 


-3.13 


-1.26 


-0.60 


-0.27 


-0.17 


50 


1.594327 


-3.38 


1.52 


-0.93 


-0.43 


-0.08 


-0.02 








1000 


2.833102 


-1.06 


0.49 


-0.28 


-0.04 















TABLE III. The double energy f%(g) = 2E*(g) of the lowest quasistationary state, as compared to numerical values f(g). 
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Ref(g) 


RefZ(g) 


-W(9) 




0.01 


0.984428 


0.984429 


0.000000 


0.000003 


0.02 


0.967451 


0.967477 


0.000001 


0.000035 


0.05 


0.900673 


0.901116 


0.006693 


0.004888 


0.1 


0.794881 


0.791404 


0.089412 


0.090883 


0.2 


0.72882 


0.728985 


0.27735 


0.281642 


0.5 


0.7477 


0.751513 


0.6100 


0.613077 


1.0 


0.8297 


0.834876 


0.9097 


0.911547 


2.0 


0.964 


0.971001 


1.260 


1.261483 


5.0 


1.23 


1.238149 


1.84 


1.832649 
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Figure Captions 



Fig. 1. The free energy of the double-well model given by the approximants f\{g) (dashed line), /2(g) (short- 
dashed line), and /2(g) (solid line). Crosses correspond to the exact values of function (61). 

Fig. 2. The real part of the double energy for the lowest quasiresonance state represented by the approximants 
Re /I (5,0) (dashed line) and Re/Kg, 0) (solid line), compared to exact values marked by diamonds. 

Fig. 3. The modulus of the imaginary part of the double energy for the lowest quasiresonance state given by the 
approximants lmf^(g, 0) (short-dashed line), Im f£(g, 0) (dashed line), and Im/g(<7,0) (solid line). Exact data are 
shown by diamonds. 

Fig. 4. The self-similar approximants /|(a, 0) (solid line) and /|(a,oo) (dashed line) for a > and the approxi- 
mants /| (a, —0) (solid line) and /|(a, —00) (dashed line) for a < 0, in the region — 1 < a < 1. 

Fig. 5. The ground-state energy for the nonlinear Schrodingcr equation (151) for the self-similar approximants 
E*(g) (solid line) and E%(g) (dashed line) and for the Thomas-Fermi approximation ExF(g) (short-dashed line). 

Fig. 6. The density (163) for the corresponding wave functions in the self-similar approximation (155) (solid line), 
Gaussian approximation (159) (short-dashed line), and Thomas-Fermi approximation (dashed line) for g = 0.2. 

Fig. 7 The self-similar approximants f£(r, 0), defined in Eq. (169), describing the structure of a vortex filament. 
The solid line is for /2*( r i0)i long-dashed line is for f£(r, 0), short-dashed line is for /|(r, 0), and the dotted line is 
for /g (r, 0) . Diamonds represent exact numerical data. 
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